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We investigate the vortex lattice solution in a (2 + 1) -dimensional holographic model of supercon- 
ductors constructed from a charged scalar condensate. The solution is obtained perturbatively near 
the second-order phase transition and is a holographic realization of the Abrikosov lattice. Below a 
critical value of magnetic field, the solution has a lower free energy than the normal state. Both the 
free energy density and the superconducting current are expressed by nonlocal functions, but they 
reduce to the expressions in the Ginzburg-Landau (GL) theory at long wavelength. As a result, a 
triangular lattice becomes the most favorable solution thermodynamically as in the GL theory of 
type II superconductors. 



PACS numbers: 11.25.Tq, 74.20.-z, 74.25.Qt 
I. INTRODUCTION 

The application of AdS/CFT (anti-dc Sit- 
ter/conformal field theory) correspondence to condensed 
matter physics has recently become one of the most 
interesting topics in string theory. (See Refs. [l|, 0] for 
reviews.) In Refs. d, 01, a holographic model of a super- 
conductor is constructed by a gravitational theory of a 
complex scalar field coupled to the Maxwell field. This 
opens up a window to study superconductors in strongly 
interacting system via AdS/CFT correspondence. 

The magnetic property of the superconductors, such as 
the Meissner effect, is a good place to start the investiga- 
tion since it plays the crucial role to distinguish between 
type I and type II superconductors Recently, many 
localized solutions under the external magnetic field B 
were constructed numerically, which include a "droplet" 
solution and a single vortex solution with integer wind- 
ing number 0, 0, i, i, [13, EH, [11 • These numerical 
solutions imply that the holographic superconductor be- 
longs to type II, in agreement with the scaling argument 
in Ref. [l|. 

There are two critical magnetic fields Bd (the lower 
critical value) and Bc2 (the upper critical value) in type 
II superconductors. At B = Bd, the external magnetic 
field begins to penetrate into the superconductor and vor- 
tices appear for B > Bd- At B = Bc2, the second-order 
phase transition occurs and the superconductivity dis- 
appears. When B approaches Bc2 from below, a vortex 
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lattice, in which a single vortex is arranged periodically, 
should appear since it is more favorable thermodynami- 
cally than a single vortex. 

In this paper, we construct the vortex lattice solu- 
tion or the Abrikosov lattice solution, which is charac- 
terized by two lattice parameters, ai and a2, perturba- 
tively near the second-order phase transition. We fol- 
low the treatment of type II superconductors based on 
Ginzburg-Landau (GL) equations. The solution includes 
a triangular lattice solution, which is known as the most 
favorable solution thermodynamically in the GL theory. 
We determine the critical value Bc2 and obtain the free 
energy parametrized by ai and a2- It is shown that the 
free energy is always smaller than that of normal state 
for any choice of two parameters. 

In the conventional Abrikosov lattice, there exists a 
circulating superconducting current around the core, and 
the current flows along the lines of constant field for the 
condensate [See Eq. (|3.33p ]. We evaluate the R-current 
to see if the current in the holographic superconductor 
has the similar property. 

It turns out that both the free energy density and the 
R-current are written by nonlocal functions of the con- 
densate unlike the GL theory. However, these expressions 
reduce to the ones in the GL theory at long wavelength. 

The plan of our paper is as follows: In Sec. [TTl we 
construct the vortex lattice solution by superposing the 
single droplet solutions found in Ref. [3]. We obtain the 
upper critical magnetic field Bc2- In Sec. lIIIl we calculate 
the free energy and the R-current in a power series ex- 
pansion in the order parameter. In Sec. llVi we take the 
long-wavelength limit. In this case, the triangular lat- 
tice solution minimizes the free energy, and the R-current 
flows along the lines of constant field for the condensate. 
Section |V] is devoted to conclusion and discussion. 
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II. PHASE DIAGRAM AND SOLUTIONS 



Eqs. (|2.4p become 



We consider a (2+l)-dimensional holographic super- 
conductor described by a dual gravitational theory in 
four dimensions (AdS^) coupled to a charged complex 
scalar field ^' and a Maxwell field To investigate 

the superconducting phase near the second-order phase 
transition, we consider the equations of motion close to 
the phase transition, as considered in Refs. [1, (3]. For 
simplicity, we take a probe limit where the backreaction 
of the matter field onto the geometry can be ignored [1] . 

In this section, we first obtain a "droplet" solution sim- 
ilar to the one obtained in Ref. Q by solving the equa- 
tions of motion at the leading order. We also obtain 
the upper critical magnetic field Bc2 as a function of T. 
Then, we construct the vortex lattice solution by super- 
posing the droplet solutions. 



A. "Droplet" solution 

The background metric is given by AdS'4-Schwarzschild 
black hole with metric 



d h d A? m^L'^ 



du V? du a^h 



L2a2 



u 



(2.5a) 



a^h^ + A]At = ^^At\^f, (2.5b) 



(2.5c) 
(2.5d) 



where i,j, k = x,y and A = + dy. 

We solve these equations under the following boundary 
conditions: 



• For the scalar field, we shall confine our interest to 
m?L?' — —2 case. Then, the asymptotic behavior 
of is 



-{-h{u)dt' + dx' +dy') + 



h{u) = l-u^ 



a{T) 



u'^h[u) 



, (2.1a) 
(2.1b) 



where L, Rq, and T are AdS radius, horizon radius, and 
the Hawking temperature, respectively. We take the co- 
ordinate u such that the horizon is located at u = 1. The 
action of the matter system S = {L? /2K\e'^)S is written 
by 



S 



d^x 



M 



p-2 

-— - -m^|*|^ ) , (2.2) 



* ~ CiU -I- C2-u^ (u->0), 

and both modes are normalizable (isj . In this case, 
both ci and C2 can be interpreted as expectation 
values of the dual operators with A = 1 and A = 2, 
respectively Q . For simplicity, we consider only the 
case where the faster falloff is dual to the expecta- 
tion value, i.e., ci = 0. We also impose that is 
regular at the horizon u = 1. 



• The asymptotic values of the gauge field give the 
chemical potential /x and the external magnetic 
field B: 



where m and e are the mass and charge of the scalar field 
^, respectively, and 



F^, ^ d,,A. - d,A^. (2.3) 



The probe limit is realized by taking the limit e ^ 00, 
keeping and ^' fixed. The equations of motion are 
given by 

- m^* = 0, (2.4a) 
V,F^" = := ^[(^?^.*)% - *^(^?^*)]. (2.4b) 

Hereafter, we choose a gauge = 0. We will con- 
sider stationary solutions since our interest is in ther- 
modynamics of the holographic superconductor. Then, 



II = At{x,u ^ 0), B{x) ^ F^y{x,u^{)). 

The boundary condition at the horizon is given by 
requiring that Af^^dx^ has a finite norm there, i.e., 
Aiix, u = 1) is regular and At{x, u = 1) = 0. 

We will only change the external magnetic field B per- 
pendicular to the AdS boundary, keeping the tempera- 
ture T and and the chemical potential /.i on the boundary 
theory fixed. In this case, the scalar field 4* begins to con- 
densate below a critical value of the magnetic field i?c2, 
while the condensation does not occur above the critical 
value. 

Defining the deviation parameter e from the critical 
magnetic field i?c2 (> 0) as e := (i?c2 — B)/Bc2, we can 
expand the scalar field ^P, the gauge field, and the current 
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FIG. 1: (color online). _Bc2(+) as a function of T/Tc- Dashed 
line: B j = (1/8) (T/r,)^ (See Sec. HV)l . 

as a series in e ^: 

^{x,u) = 6^/^1(3;,") + e3^V2(ic, «) + ••• , (2.6a) 

A^,{x,u)^Af{x,u) + eA^^\x,u) + -- - , (2.6b) 

j^{x,u) = ejl'^x,u) + e^jj^Hx,u) + --- , (2.6c) 

where x = {x,y). 

The zeroth order solution generating the critical ho- 
mogeneous magnetic field Bc2 (> 0) and the chemical 
potential /i are given by 

^i°)=Ml-«), v4f=0, 4°)=i3e2X. (2.7) 

Substituting Eq. (|2.7p into Eq. (|2.5ap and taking an 
ansatz ipiix^u) = e'™^(a;, u;p)/L for a constant p, the 
equation of motion for ipi reduces to 

2 d (h{u) d \ {A^t\-^)f m^L^ 



du \ du 



(t>{x,u;p). 



(j){x,u;p) 



(2.8) 



Following the ansatz in Ref. Q, we separate the vari- 
able (f) as 0„(x, u;p) = p„ {u)^„ {x; p) / L with a separation 
constant A„. The equations for p„ and 7„ are divided into 
the following equations: 

^ ^ 7„(x;p) 4r7n(a;;p), 



4 
2h 
u 
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1^ /I \2 . 2 Bc2\ 



(2.9a) 



(2.9b) 



where X := y''2Bc2{x ~ p/Bc2) and g := /i/a is a dimen- 
sionless quantity. 

The regular and bounded solution of Eq. (|2.9ap satis- 
fying lim|^|^oc l7ri| < 00 is given by Hermite hmctions 
Hn as 



7„(x;p) =e-^'/'*i/„(X), 
and the corresponding eigenvalue An is 
A„ = 271 + 1, 



(2.10) 



(2.11) 



for a non-negative integer n. 

This solution corresponds to a "droplet" solution ob- 
tained in a earlier work 0] in the sense that these so- 
lutions fall off rapidly at large A single (localized) 
droplet solution is easily obtained by generalizing the so- 
lution with circular symmetry, as shown in Refs. 

Now, let us consider the phase diagram. Equa- 
tion (j2.9bp is characterized by two parameters T//x oc 1/q 
and Bc2^n/tJ'^- Also, recall that we impose boundary 
conditions both at asymptotic infinity and at the hori- 
zon. Such a problem has a nontrivial solution only when 
there is a relation between these two set of parameters. 
The solution corresponds to the case where the 'J''-"-' = 
state becomes marginally stable. 

If one has a second-order phase transition from the ^' = 
state to the "hairy" black hole, one should encounter a 
marginally stable state at the point of transition. So, the 
solution of Eq. (j2.9bp gives the candidate for the phase 
transition point. 

It is clear which side of the phase transition line is 
the superconducting phase. Suppose that one lowers the 
magnetic field with a fixed T/fi. The right-hand side 
of Eq. (j2.9bp suggests that the magnetic field increases 
(effective mass)^ and tends to stabilize the state. [On the 
other hand, the electric field decreases (effective mass)^ 
and tends to destabilize the 5*^°) = state.] Thus, the 
xpC) — state should be stable under a large enough 
magnetic field, and the state becomes unstable as one 
lowers the magnetic field. Then, it is likely that the up- 
per critical magnetic field is given by the value of the 
magnetic field when one first encounters a marginally sta- 
ble state (as one lowers the magnetic field). Note that 
a marginally stable solution is parametrized by BXn/p^ 
for a fixed T/ fi. Thus, one has the largest magnetic field 
when A„ takes the minimum, namely the n — solution. 

Consequently, the upper critical magnetic field Bc2 is 
given by the largest B/n^ for the solution of Eq. (|2.9b|) . 
Figure [T] shows Bc2 obtained in this manner in the 
{T/Tc, B / fi^) phase diagram.^ Here, Tc is the critical 
temperature when there is no magnetic field, and it is 



^ This type of expansion has been applied to the SU(2) model of 
the holographic superconductor and the speed of second sound 
was analytically derived [Tsll . 



^ Obviously, the above argument gives only the necessary condition 
for the phase transition and does not give the sufficient condition. 
The sufficient condition is given by showing that there exists a 
condensate solution for B < Bc2 and that its free energy is lower 
than the 'I' = solution. This is shown perturbatively below. 
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determined in the combination Tc//i = 3/(47r(7c) with 



B. Vortex lattice solution 

Let us clarify how the vortex lattice constructed in this 
subsection differs from the droplet in the previous sub- 
section. A characteristic feature of a vortex is that it 
has a zero and has a winding number around the zero. 
Thus, the droplet which is nonvanishing everywhere can- 
not have such a winding number. However, a superposi- 
tion of droplets can have zeros and winding numbers as 
we will see in a moment. 

As we have seen in the last subsection, it is enough to 
consider only the n = solution near Bc2- 



lo{x;p) = e 



exp 



1 f X 

o pro 

2 V^o 



(2.12) 



where tq := l/\/Bc2- As A,i is independent of p, a linear 
superposition of the solutions e*P^p„(M)7„(x;p) with dif- 
ferent p is also a solution of the equation of motion for 
^ at 0(e^/^). To obtain the vortex lattice solution from 
the single droplet solution (|2.12p . consider the following 
superposition: 



L 



J2 cie'P'y-fo{x;pi), 



ci 



exp 



/ — — oo 



-I' 



at 



Pi 



2ttI 



(2.13a) 
(2.13b) 



for arbitrary parameters ai and 02. In terms of the ellip- 
tic theta function i)^ defined by 




FIG. 2: The vortex lattice structure for the triangular lattice 
in the (a;, y)-plane. The vertical line represents a = |7i|^ and 
vortex cores are located at |7l| = 0. 



so Tpi also has a pseudo-periodicity 

i)i{x, y, u) = 'ipi{x,y + aiTQ, u), 



(2.18a) 



27rro a2ro 

Vi[x^ ,y-\ ' " 



ai 



ai 



exp 



"27ri 






ai 




2a J _ 



Mx, y, u). (2.18b) 



Thus, <7{x) :— |7L(a;)p represents a lattice in which the 
fundamental region Vq is spanned by two vectors 61 = 
o-irody and ^2 = 2'Kr^laidx -\- a^T^jaxdy and the area 
is given by 27rrQ. This is the well-known result, where 
the magnetic flux penetrating the unit cell is given by 
Bc2 X (Area) = 27r. This shows the quantization of the 
magnetic flux penetrating a holographic vortex. 
Second, -dz vanishes at 



1,. 



(2.19) 



(9 



e""), (2.14) 



/— — 00 



the summation over I in Eq. (|2.13ap is expressed by 



(2.15) 



/— — 00 



where 



-ix + y 
aifo 



2i:i — a2 



(2.16) 



The solution (|2.13p or (|2.15p represents a vortex lattice. 
The elliptic theta function i?3 has two properties which 
are useful to see the vortex lattice structure. First, 
has a pseudo-periodicity 

dz{v + \,T)^'dz{v,T), (2.17a) 
dz{v + T, r) = e-2-*(''+"/2) dz{v, r), (2.17b) 



for any integers m, n. Since the expectation value of the 
operator O dual to ^' is proportional to 71,(3;), the con- 
densation (O) has Qu zero at Xjyi n. Also, it is easily shown 
that the phase of (O) oc rotates by 27r around each 

Xm,n [El- Thus, the cores of vortices are located at x,n^n- 
The triangular lattice, where three adjoining vortices 
Xm,n form an equilateral triangle, is given by the follow- 
ing parameters: 



02 
ai 



Ol 

2 



3-i/S 



(2.20) 



Figure shows the conflguration of <j{x) = \'^l{x)\'^ in 
the (cc, y)-plane for the triangular lattice. It is well-known 
that this conflguration minimizes the free energy in the 
GL theory [|. 

Obviously, the linear analysis presented here alone can- 
not lift the degeneracy of the solution and determine that 
the triangular solution is the correct one. To determine 
the correct conflguration, one needs to compute the free 
energy and include nonlinear effects for the holographic 
superconductor as well, which is a goal of later sections. 
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III. FREE ENERGY AND R-CURRENT 

111 this section, we calculate the free energy of the vor- 
tex lattice solution derived in the previous section. For 
simplicity, we consider the system in which the region V 
for scalar field condensation is compact and very large 
compared with the unit cell of the lattice. So, we con- 
sider the system in which |^| is zero outside the region V. 
Since we fix both of the chemical potential ^ = At{u ^ 0) 
and the temperature T under the variation of the mag- 
netic field, A'i'\u = 0) = = 1, 2, • • • ) throughout this 
article. 



Preliminaries 



Employing the equations of motion (|2.4ap . the on-shell 
action from Eq. p.2[l becomes 



S --i 



-gF^ 



M 



dM 



(3.1) 



where dS^ is the integral measure normal to the space- 
time boundary dAi. 

The boundary d^A in Eq. (j3.ip consists of 



• two spacelike surfaces (the past and future sur- 
faces), 

• the horizon, 

• the AdS boundary dM.i, 

• the boundary in the [x, y)-coordinates dM2 {x = 
const, and y = const, surfaces). 

We consider a stationary problem, so the contribution 
from two spacelike surfaces cancels each other. There is 
no contribution from the horizon as well from the reg- 
ularity condition at the horizon and from the "rcdshift 
factor" hill ^ 1) = 0. Thus, we need to consider the 
remaining two surfaces dMi and dM2 to evaluate the 
on-shell action. 

Furthermore, one can ignore the second term in 
Eq. ()3.1|) because the scalar field 4" has a compact sup- 
port in the (x, i/)-coordinates and it satisfies the bound- 
ary condition ci = at the AdS boundary. Therefore, 
the on-shell action reduces to^ 



S --1 



d^X 



M 



Let us expand the on-shell action 



as 



(3.2) 



(3.3) 



^ In principle, one also should take the counterterm action Set for 
the scalar field ^ into account, but one can easily see that 5ct 
makes no contribution under the boundary condition at the AdS 
boundary. 



In terms of gauge field strength pjill defined by 

d^A^ - d,Af (* = 0, 1, 2, • • • ), (3.4) 
the coefficient S^^^ becomes 



-5(1) = 



d^x 



M 



'_Ji pM'' _ 

2 ^(0) ^t^i' ^ 



dM 



dii,Fi;^A\^\ 



(3.5) 



where we used the Maxwell equations PjS) ~ 0. 



Similarly, S^"^^ becomes 



5(2) = - 



dM 



^(0) 



X v-g J, 



M 



_ fA"' 4(1; 



where we used the Maxwell equations V'^pj^^ = and 

"^"F^D = jI^^ The last term of Eq. dHU) vanishes from 
the "orthogonality condition" given in Eq. (|A3[) . Thus, 
the on-shell action up to second order in e is given by 
boundary terms in Eqs. p.5p and (j3.6p . 

As mentioned above, wc need to consider two surfaces 
dA4i and 9A^2-'* However, our interest is in the config- 
uration of ^ which gives the lowest free energy, namely 
the ^-dependence of the free energy. Thus, it is enough 
to extract only the ^'-dependence of the on-shell action. 
When ^' has a compact support in the {x, ?/)-coordinates, 
dM2 gives the terms which does not depend on ^f, so 
they are irrelevant.^ Consequently, the boundary terms 
we arc interested in entirely come from dMi and the 
on-shell action up to second order in e is given by 



(3.6) 



- — 



dMi 



-^(0) ^1^+2 (1) " 



d^x 5''F},,> A 



(1) 4(1) 



+ 0(e3). 



O(e^) 
(3.7) 



u=0 



Note that one can ignore the first term in the first line 
since it does not depend on 5". Also, we used ^j^'' {x, 0) = 
0. 

According to the AdS / CFT dictionary, the expectation 
value of the R-currcnt in the boundary theory, ( J^{x) ) 
is calculated by 



aFui_,{t,x,u 0). 



(3.8) 



* The surface dA42 is often negligible in thermodynamic limit, 
but this is not the case in the presence of a magnetic field [See 
Eq. ITTt ], 

^ The field has a compact support in the (x, j/)-coordinatcs. But 
the gauge field has a noncompact support which depends on 
ipi [See Eq. Il3.15l l1 . Thus, one would have ^t-dependent terms 
from dM2 through the gauge field. However, the i/)i-dcpondent 
part of A/^ decays exponentially as seen in Sec. IIV Al 
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Since f!? = 0, 



aeFi]\x,u^0) + O{e^). (3.9) 



One can rewrite the on-shcU action (|3.7p by the R- 
current expectation value (j3.9p as 



(1) 



+ 0(e3). (3.10) 



'U=0 



Since the free energy Q is related to the Lorentzian on- 
shell action as $7 = —Sos/ J dt, 



(1) 



+ 0(e3). (3.11) 



11=0 



Note that these expressions ignore the terms which do 
not depend on 5". 

Both for the free energy and for the R-current, the 
problem amounts to evaluate ( ) cx F^^\ so we obtain 
this quantity in the next section. 



B. First order solution of Au 



The Maxwell equations at first order are given by 



regular at the horizon and eFxl^ = B — Bc2 at the AdS 
boundary, namely, 2d[xA^y^ {u = 0) = — -Bc2- 

The solutions are formally written by introducing 
Green functions which satisfy the Dirichlet condition at 
the AdS boundary: 



A\ 



/o u'2/l(u')' 

X j dx'Gt{u]u\x-x')a{x'), (3.15a) 



A^>=a,{x)~a'e,= du' ^ 

X j dx'GB{u;u'\x-x')djG{x'). (3.15b) 

Here, ai{x) is a homogeneous solution of Eq. (|3.13bp sat- 
isfying 2d[xO'y] = —Bc2 at the AdS boundary, and is 
independent of w. Gt and Gb are Green functions of 
Eqs. (|3.13ap and (j3.13bp . respectively: 

h — + A] Gt(u; u'lx) = -h{u) 5(u - u')S(x), 
ow' J 

(3.16a) 

Gt{u = 0; u'\x) = Gt{u = 1; u'\x) = 0, (3.16b) 



92 



(1) _ 2L2q,2 



a'hj-^ + AjAY' = —^An^^\', (3.12a) 



1) 



0. 



(3.12b) 
(3.12c) 



We choose the gauge A„ = 0, but there is a residual gauge 
transformation Ai — > Ai — diA{x). From Eq. (|3.12c[) . one 
can set S'3^^Af^ = using the residual transformation. 
{ Jfi) oc Fu^ is of course invariant under the transfor- 
mation. In this gauge 8^^ diA^'p = 0, Eq. (|3.12[) becomes 



a^h 1^ + a) A^P ^ ^^aTo[x\ (3.13a) 



d d 



\ 2 2 

^a'^h-^ + A] = ^ e,^d,aix). (3.13b) 



Here, is the antisymmetric symbol ei2 = —£21 = 1 



and e;^ = dj. We also used 



■ d) ^ Pl{u)d(y{x) (1) _ pI{u) da{x) 



from Eqs. ((Ob)) and ([^1^ . 

The boundary conditions for Eq. (|3.13ap are = 

4^^(u = 1) = = 0)- For Eq. (|3.13bp . we im- 

pose the boundary conditions such that the solution is 



du V d 



-A Gb{u-u'\x) ^ ^5{u-u')5{x), 

(3.17a) 

Gi3(M = 0; u'\x) = lim /i(u)(9„G'b(w; u» 0. (3.17b) 



The boundary condition Eq. (j3.17b[) at it = 1 follows 
from the regularity at the horizon. Similar expressions 
hold for u ^ u' from G{u] u') ~ G{u'; u). 



C. Free energy and R-current 

Using the Green function Gb, one can express the R- 
current expectation value as 



Ji{x) ) = -eidjE{x), 



(3.18) 



where 



E{x) 



dx'a{x') 



X du C du'^^^GB{u;u'\x-x') 
Jo w 



(3.19) 



Here, we used OGb/Ox'^ = -dGB/dx\ Equation (jXTBl) 
implies a circulating R-current which is dictated sim- 
ply from the current conservation. It is a different issue 
though whether the current actually flows along the lines 
of constant field for the condensate like the GL theory. 
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The free energy can be expressed by the Green function 
similarly. The free energy is given by 



using Eq. (PTTjl and 2d[^A'-y^^\i 
Eqs. (jXTT)) . Gb satisfies 



0) 



-Bc2- From 



2 d f d 



du V du 



/ da; Gb{u; u'\x) ~ 

JR2 



= —d{u — u'), 



and 



/ dxGBiu;u'\x) 



min(«,«') 



Consequently, one obtains 



Cvo\{V)a 



where 



(3.20) 



(3.21) 



(3.22) 



and / indicates the average of / over V in the {x,y)- 
plane: 



1 



vol(V^) 



dx f{x). 



(3.23) 



From Eq. (|3.2ip . < 0, which suggests that the free 
energy of the vortex lattice state is smaller than the one 
of the normal state (4" = 0). 

Equations ((3T8)) . ((3J9)) and ((3^ arc not our final 
results. Because po and 7^ arc the solutions of the linear 
equations (|2.9bp and (|2.9ap . their normalizations have 
not been determined yet. Below wc eliminate this am- 
biguity, so ( ) and arc characterized only by two 
lattice parameters ai and 02- 

The ambiguity in the normalizations come from the 
linear equations, and it is resolved only after one con- 
siders nonlinearity. So, let us consider the equation of 
motion for -02. As shown in Appendix A, if ip2 obeys the 
boundary condition ci = and the regularity condition 
at the horizon, we obtain the "orthogonality condition" 
(|A3|) which relates po to the gauge field A^^^ . 

Using Eqs. (|3.14p and integration by parts, the "or- 
thogonality condition" (|A3p can be written as 



du 



/ dxF^^^a{x) 



= 2 



duJ^A^\u) 
u^h[u) 



dx (^{x). (3.24) 



From Eq. (|3.15bp . 

i^ij) = -Sc2 + du'P^A[GB{u; u') * a]{x), 

Jo " 

(3.25) 

where " * " is convolution in the (x, ?/)-plane: 

[f*9]ix):= j dx'f{x-x')g{x'). (3.26) 



Substituting the solutions p.lSap and (|3.25|) into 
Eq. (1X211), wc obtain 

B.2C C d^PM C d^' pMIiM^ (3.27) 



O- JO 

where we define 







w 



a 



24"^ (w) a{Gt(u;u')^a\ 2A'^\v:) 



1 : 



h{u) 



h{u') 



a A Gb(u; u') * a 



(ay 



(3.28) 



which is independent of the normalization of 7^ . 
Then, the free energy (|3.2ip is expressed by 



L2 625^2 a vol(y) 



T 



(3.29) 



where 
r := 



in) f\^^Pl^j^,y) 



du 



■I{u,u). (3.30) 



/o Jo " 

Here, po is the solution of Eq. (|2.9bp normalized by 



du ^ = 1. 



(3.31) 



Since the functions T and F do not depend on the nor- 
malization of Po and 7l, the free energy (|3.29p also does 
not depend on the normalization cither. 

Similarly, the potential S (|3.19p which gives the R- 
current expectation value does not have the ambiguity 
in normalization. This can be seen by expressing the 
potential as 



E{x)^eB, 



L-'a^ 1 



2^262 r 



du 



\.^ rpl{n') [GB{uW)*(r]{a 



(3.32) 



tt=0 



This is our main result: the free energy (|3.29p and the 
expectation value of the R-current (|3.18p are expressed 
by F (IX^ and S 
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As seen in these expressions, the free energy density 
and the R-current are not only determined by the com- 
plex scalar field a — at x but also by <t in the entire 
region around x. On the other hand, in the GL theory, 
the free energy density is a local function of the order 
parameter "^qi^^x). The superconducting current is also 
expressed locally by 



J^'^{x) = -e^d,\^GL{x) 



(3.33) 



This implies that the superconducting current flows along 
the lines of IvPglP = constant Q. These differences are 
natural since we have not taken a long-wavelength limit 
in evaluating the free energy and the R-current unlike 
the GL theory. By taking a long- wavelength limit, our 
result should reduce to the GL theory. As a result, the 
triangular lattice solution, which is the most favorable 
solution in the GL theory, should also become the most 
favorable one in the holographic superconductor. One 
should reproduce Eq. (|3.33[) as well. We will demonstrate 
this in the next section. 



IV. THE LONG- WAVELENGTH LIMIT AND 
THE TRIANGULAR LATTICE 

A. A representation of Green functions 



The expression of F in Eq. p.30p is nonlocal both in the 
AdS radius direction and in the (x, ?/)-directions, which 
makes the boundary interpretation rather unclear. The 
nonlocalities come from Green functions Gb and Gt, so it 
is useful to expand the M-dependence of the Green func- 
tions by a complete set of orthonormal functions. 

First, let us consider eigenfunctions xa(w) satisfying 

CbXx{u) = AxaH, Cb := -a^-^ ( ^T" ) ' (^.la) 

au \ du I 



Xa(w = 0)= \iuih{u)^^{u)^Q. 



(4.1b) 



The operator Cb becomes an Hermitian operator for the 
inner product defined by 



(j) \ iP)b ■= du (j)^u)ip{u) 
Jo 



(4.2) 



for any solutions (f) and ip of Eq. (14. ip . One can easily 
show A > 0. 

The normalized eigenfunctions {xa} with respect to 
the inner product (j4.2p form a complete orthonormal set 

(xa|xa'>s = (^AA', ^xx{u)xl{u') ^ 6{u- u), (4.3) 



and the Green function Gs in Eq. (|3.17p is represented 
by {xa} as 



Gb{u; u'\x) = J2 Xx{u)xliu') G2{x- A). 



(4.4) 



Here, G2{x\xv?) is the Green function defined on {x,y)- 
plane as 

(A - m'^)G2{x; m^) ^ -6{x). (4.5) 

The solution of Eq. (j4.5p satisfying hm|a,|^oo IG2I < 00 
is given by the modified Bessel function: 

G2{x;m^) = ^Ko{m\x\) (4.6) 

for any real positive value of m^. 

Similarly, Gt in Eq. p.l6p is constructed as 

Gt{u;u'\x) = Y,U^'Kli^')G2{x;v), (4.7) 

r;>0 

where {^j;} is the complete orthonormal eigensystem of 
the equation 



du? 



C^(« = 0)=e^(u=l) = 0, 



(4.8a) 
(4.8b) 



with 



The inner product ( 



)t is defined by 
du 



K^) 



4>\u) ^(u) 



(4.9) 



(4.10) 



for any solutions </> and ip of Eq. (|4.8p . The operator Ct is 
Hermitian with respect to the inner product (|4.10p . and 
?7 > 0. 

Substituting Eqs. dM]) and gJl) into Eq. ([3^ . we 
obtain another expression of F as 



F = ^P(A)Ci(A)+^0(,7)Co(?7) 

A V 

where we introduce Cri('Ti^) 
C«(m2) 



2, crA"[G2(m2) *cr 



and the functions P, Q defined by 



P(A) := 



Qiv) ■■= 



(4.11) 

(4.12) 

(4.13) 
(4.14) 



Similarly, ^ (|3.32p is represented as 



E{x)^eB,2 



Exa(O) 



2nle^ F 



A>0 



(4.15) 
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The eigenvalues A and rj can be obtained numerically 
by solving the two differential equations (|4.ip and (|4.8p . 
Both are positive, and the minimum values of A and ?/ 
are given by A ~ 2.22a^ and 77 ~ 7.41a^, respectively. 

The holographic superconductor is constructed in the 
gravity theory with one extra dimension which is ex- 
tended perpendicular to the (2 + l)-dimensional space- 
time of the boundary theory. So, from the (2 -|- 1)- 
dimensional point of view, the bulk gauge fields with 
a wide variety of mass \/A and appear as in Eqs. (|4.ip 
and g^. 



B. The long- wavelength limit 

We find that both the free energy density and the R- 
current take nonlocal forms whereas they take local forms 
in the GL theory. This is because the AdS/CFT results 
correspond to the results to all orders in effective the- 
ory expansion. The GL theory takes only first few terms 
in the effective theory expansion, so one needs to take a 
long- wavelength limit in the AdS/CFT results to com- 
pare with the GL theory. Taking such a limit is common 
in hydrodynamic studies based on the AdS/CFT duality 
(See, e.g., Refs. [13, H H-) 

In this subsection, we demonstrate that our results in- 
deed reduce to the GL theory ones in the long-wavelength 
limit. As a result, the triangular lattice solution becomes 
the most favorable one like the GL theory. 

We have obtained nonlocal expressions, but there is in 
fact an analogous situation in a superconductor. A su- 
perconductor has small length scales, the Pippard/BCS 
coherence length and the mean-free path. In the pres- 
ence of these small scales, the electromagnetic response 
in general takes a nonlocal form. The local form such 
as the London equation is the limit where these length 
scales are negligible Hf . 

In our problem, the nonlocalitics come from the con- 
volution in Eqs. (|iT^ and ((iT5)l . i.e., 

[G2(A) *cr](a;) = [ dx'G2ix ^ x'; X) a{x') (4.16) 



where we used 



and a similar expression for G2{ti). The Green's functions 
G2 have the natural length scales 1 / \/A and 1 / They 
are the small length scales in our problem. Their micro- 
scopic interpretation is unclear but they are 0{T~^). On 
the other hand, the natural length scale of the vortex lat- 
tice is the size of the fundamental region parametrized by 
tq. This is the length scale of the condensate a = |7lP- 
When To 3> 1/ \/A, 1 / y^, G2 quickly decays compared 
with cr, and the convolution reduces to a local form. This 
is what we meant by the "long-wavclcngth limit." 

By replacing a{x' ) by <7{x) in the convolution, we ob- 
tain 



'G2[X)*(j][x) = 



(7{x) 



1 + 



1 



(4.17) 



/ dx G2{x; A) = Y- 

JR2 A 



Subleading terms in Eq. (|4.17p can be estimated from a 
series expansion a{x') = a{x) -I- • • • . 

Substituting Eq. (|iTT7)) into Eq. (|iT^ . one obtains the 
coefficient F (|4.1ip as 



r ^ 



1 1 

rgry' r^A 



1 + 



where 



?;>0 



Qiv) 

V 



(4.18) 



(4.19) 



r > from its definition ()3.30p . [The constant Ci is 
indeed positive from Eq. (|4.14p and the positivity of rj.] 
One has the lowest free energy D, when F takes its mini- 
mum. Therefore, the thermodynamically realized config- 
uration is given when 



(a)2 



(4.20) 



takes the minimum value. This is the same condition 
as the one for the Abrikosov lattice in standard type II 
superconductors As is well-known, the minimum is 
P ~ 1.16, which is given by the triangular lattice (|2.20p . 
Similarly, the R-current (|4.15p becomes 



E{x) = eBc2 



L'^a 1 a{x) 



(4.21) 



using Eq. (|4.17p . Thus, like the GL theory, the circulat- 
ing R-current flows along the lines of constant field for 
the condensate. 

Let us consider the validity of our approximation. One 
can achieve the long-wavelength limit by taking a small 
magnetic field Bc2 since rg = As shown in 

Fig. [TJ one can take an arbitrary small value of Bc2 by 
choosing the temperature T suitably. Using q = /i/a, A ~ 
2a^, and Qc ^ 4, the condition Bc2 <C A, 77 gives 



2 2 

M 9c 



(4.22) 



This condition is implemented in Fig. [T] Our approxi- 
mation in this subsection becomes good if the system is 
located well below the dashed curve in the phase diagram. 

Incidentally, one would incorrectly conclude that the 
expanded solution of Eqs. (|2.6p exists even for B > Bc2 
by replacing the deviation parameter with e = [B — 
Bc2)/Bc2 (> 0). Suppose that this is the case. Then, 
this changes the sign in the first term of FiJ' in the so- 
lution (j3.25p . As a result, one has an extra minus sign in 
the following equations: the left-hand side of Eq. (|3.27p 
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and the right-hand side of the second line of Eq. (|3.30p . 
But this contradicts with the result in Eq. (|4.18p . There- 
fore, the e-expansion of ^ does not exist for B > Bc2 at 
least when rg ^ 1/ V\, 1 / y^. This agrees with the ar- 
gument in Sec. II that the superconducting phase ceases 
to exist at B = Bc2 and the normal phase appears for 
B > B,2. 



V. CONCLUSIONS AND DISCUSSION 

We have investigated the vortex lattice solution in the 
holographic superconductor. One would interpret the 
holographic superconductor as a superfluid, and the com- 
putations described here can equally apply to this case 
as well in a slight modification. In this case, the rotation 
of the superfluid is analogous to the magnetic field. 

One main difference between the conventional super- 
conductors in the GL theory and the holographic super- 
conductor is that both the free energy density and the 
R-current are expressed by nonlocal quantities, such as 
the two-point function G2 in Eq. (|4.5p . This is because 
our results in Sec. IIIII is beyond the applicability of the 
GL theory. In fact, we show in Sec. [IV] that our results 
reduce to the results of the GL theory by taking a long- 
wavelength limit. Then, one should be able to estimate 
the corrections to the GL theory by inspecting our ex- 
pressions in Sec. IIIII 

Another interesting direction to pursue is the investi- 
gation of the kinetics of the vortex lattice solution. For 
example, it is well-known that vortices move at a con- 
stant velocity along the direction perpendicular to the 
magnetic field by the Lorentz force and the electric resis- 
tance appears if we add an appropriate external electric 
field perpendicular to the magnetic field. This phenom- 
ena is important for pinning of superconductors with ir- 
regularities 0. We will discuss this in more detail in 
Ref. [l^]. It is also interesting to consider the dynamic 
critical phenomena of the vortex lattice solution, as in- 
vestigated in Refs. [HI, HI]. 
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APPENDIX A: "ORTHOGONALITY 
CONDITION" 

The equations for ipi and 'ip2 are given by 



Z?^o) - "I'j V'l = , i^Dl,^-m'jij2^J , (Al) 
J ^ {Z?f (a^^) ^1) + , (A2) 



where D'^^^ 



Recall that the field ^' has a compact support in the 
(x, ?/)-coordinates and satisfies the regularity condition 
at the horizon and ci = at the AdS boundary. Then, 
using Eq. (|Aip and integration by parts, one obtains the 
"orthogonality condition:" 



OM I J 

/ . (A3) 
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